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ABSTRACT 

Multi-frequency radio observations of the radio bridge of powerful classical double 
radio sources can be used to determine: the beam power of the jets emanating from 
the AGN; the total time the source will actively produce jets that power large-scale 
radio emission; the thermal pressure of the medium in the vicinity of the radio source; 



> 

and the total mass, including dark matter, of the galaxy or cluster of galaxies traced 

in' 

t^. ■ allow a determination of each of these quantities using radio observations are presented 



by the ambient gas that surrounds the radio source. The theoretical constructs that 
allow a determination of each of these quantities using radio observations are presentee 
and discussed. Empirical determinations of each of these quantities are obtained and 



o 
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Q\ . analyzed. 

c| . A sample of 14 radio galaxies and 8 radio loud quasars with redshifts between 

zero and two for which there is enough radio information to be able to determine 
the physical parameters listed above was studied in detail. The beam power was 
determined for each beam of each AGN (so there are two numbers for each source), 
and these AGN are highly symmetric in terms of beam powers. Typically, the beam 
powers are about 10 erg s . No strong correlation is seen between the beam power 
and the core - hot spot separation, which suggests that the beam power is roughly 
constant over the lifetime of a source. The beam power increases with redshift, which 
is significant after excluding correlations between the radio power and redshift. The 
relationship between the beam power and the radio power is not well constrained by 
the current data. 

The characteristic or total time the AGN will actively produce a collimated outflow 
is determined. Typical total lifetimes are ~ 10 7 years. Typical source lifetimes decrease 
with redshift, and this relation can explain the decrease in the average source size 
(hot spot - hot spot separation) with redshift. Thus, high-redshift sources are smaller 
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because they have shorter lifetimes, since the rate of growth of sources is found to 
increase with redshift. 

A new method of estimating the thermal pressure of the ambient gas in the 
vicinity of a powerful classical double radio source is presented. This new estimate 
is independent of synchrotron and inverse Compton aging arguments, and depends 
only upon the properties of the radio lobe and the shape of the radio bridge, which 
is used to determine the Mach number of lobe advance. A detailed radio map of the 
radio bridge at a single radio frequency can be used to estimate the thermal pressure 
of the ambient gas. Thermal pressures on the order of lCP 10 dyne cm -2 , typical of 
gas in low-redshift clusters of galaxies, are found for the environments of the sources 
studied here. It is shown that appreciable amounts of cosmic microwave background 
diminution are expected from many of these clusters, which could be observable at 
high frequency where the emission from the radio sources is weak. 

The total gravitational mass of the host cluster of galaxies is estimated using the 
composite pressure profile and the equation of hydrostatic equilibrium for cluster 
gas. Masses and total density profiles very similar to those of low-redshift clusters 
of galaxies are obtained. Thus, some clusters of galaxies exist at redshifts of two or 
so. The redshift evolution of the cluster mass is not well determined at present. The 
current data do not indicate any negative evolution of the cluster mass, contrary to 
what is expected in a high density universe. 

Subject headings: galaxies: active — quasars: general — radio continuum: galaxies 

1. Introduction 

Classical double radio sources, also known as FRII radio sources (Fanaroff k, Riley, 1974), 
have been extensively studied. Numerous surveys of radio sources as well as detailed studies of 
individual sources have provided a wealth of radio data on many sources. Theoretical models for 
the sources developed and discussed by various authors (e.g. Blandford & Rees, 1974; Scheuer, 
1974, 1982; Begelman, Blandford, & Rees, 1984; Begelman & Cioffi, 1989; Daly 1990) lend much 
understanding to the physics of the sources. It is believed that FRII sources are powered by highly 
collimated outflows from an active galactic nucleus (AGN), and that the radio emission from these 
sources is the result of an interaction between the beam, or jet, and the ambient medium. As a 
result, careful study of FRII sources can yield very useful information about the FRII sources, 
their environments, and their central engines. 

Daly (1994, 1995) presents a model for powerful extended FRII sources and showed how key 
parameters of an FRII source and its environment, such as the beam power and the ambient gas 
density, can be estimated using multi-frequency radio data. Results for a sample of powerful 3CR 
radio sources are also presented by the same author (Daly 1994, 1995). A larger sample was later 
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compiled by Daly and collaborators and is studied in detail. Many results from these studies have 
been presented in a series of papers (Wellman & Daly 1996a,b; Wellman 1997; Wellman, Daly, & 
Wan 1997a, b [hereafter WDW97a,b]; Wan & Daly 1998a,b; Guerra & Daly 1996, 1998; Guerra 
1997; Daly, Guerra, k Wan 1998; Guerra, Wan, & Daly 1998; Wan 1998). For example, detailed 
results on the density and temperature of the ambient medium of the FRII sources are presented 
by WDW97a,b, and the redshift evolution of the characteristic size of an FRII source and its 
application for cosmology are presented by Guerra & Daly (1996, 1998). 

This paper presents results on several key parameters of the FRII sources and their 
environments that were not presented and discussed in previous papers. These parameters include 
the luminosity in directed kinetic energy of the jet, also known as the beam power, Lj, the 
total time the AGN will produce the collimated outflow, t*, the thermal pressure of the ambient 
medium, P t h, the effect of the hot gas on the microwave background radiation, and the total 
gravitational mass of the host cluster (we find that the radio sources lie at the centers of gas rich 
clusters of galaxies). The way that these parameters can be estimated from radio data, and the 
empirical results for the sample are presented here. 

The paper is structured as follows. A brief description of the sample and data for the sample 
are presented in §^. Discussions on the jet luminosity and timescale of jet activity are given in 
§||, the thermal pressure of the ambient medium and estimates of the effect of the hot gas on the 
microwave background radiation are presented in §|J and estimating the total gravitational mass 
of the the host cluster of the radio source is discussed in §[|. Each of these sections are divided 
into two subsections, with the first subsection presenting theory and the second one presenting 
empirical results for the sample. Further discussions of the results and a summary of the paper 
are presented in §^. 

Values of Hubble's constant H Q = 100 h km s~ 1 Mpc~ 1 and the de/acceleration parameter 
g o =0 (an open empty universe with zero cosmological constant) are used to estimate all the 
parameters in this study. Different choices of cosmological parameters, within reasonable limits, 
give results that are very similar to, and consistent with, those presented here. 



2. Sample Description and Data 

The sample used in this study has been described in detail in several previous papers (e.g. 
WDW97a,b; Guerra & Daly 1998). The readers are referred to these papers for a full description of 
the sample. A brief description of the sample is as follows. All the FRII sources in this sample are 
very powerful 3CR sources, having 178 MHz radio power greater than 3 x 10 26 h~ 2 W Hz~ sr -1 . 
They are drawn from the samples of Leahy, Muxlow, & Stephens (1989; hereafter LMS89), which 
contains sources with large angular sizes, and Liu, Pooley, & Riley (1992; hereafter LPR92), which 
contains sources with smaller angular sizes. The final sample contains 27 radio lobes from 14 FRII 
radio galaxies, and 14 lobes from 8 FRII radio loud quasars. These sources have redshifts ranging 
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from zero to 2, and projected core-hot spot separations between 25h 1 kpc and 200/i 1 kpc. Note 
that the questions of projection effects and radio power selection effects have been addressed in 
great detail by Wan k Daly (1998a,b). 

WDW97a,b used the radio maps from LMS89 and LPR92 to estimate the width of the radio 
bridge behind the radio hot spot, cll, the non-thermal pressure of the radio lobe that drives the 
forward shock front, Pl, and the rate at which the bridge is lengthening, referred to as the lobe 
propagation velocity vl- The radio information may also be used to estimate the beam power and 
total time the source will be active, as described in §3. 

In order to estimate the Mach number of lobe advance, which in turn may be combined with 
the lobe propagation velocity vl to estimate the temperature of the ambient gas, high quality 
maps which image large portions of the radio bridge are required; the radio bridge is defined as 
the radio emitting region that lies between the radio hot spot and the origin of the host galaxy or 
quasar. There are 16 bridges in the sample for which we have a detection of the Mach number 
of lobe advance, which includes 13 radio galaxy bridges and 3 radio loud quasar bridges. Thus, 
the ambient gas temperature and pressure may be estimated for these sources. Lower bounds on 
the Mach number, and hence upper bounds on the ambient gas temperature and pressure, are 
available for the 14 other bridges, including 8 galaxy bridges and 6 quasar bridges, as described 
in WDW97a,b. The maps of the sources do not have sufficient dynamic range to allow detailed 
analyses of their bridge structure, and thus do not have estimates on the Mach number of lobe 
advance, nor on the ambient gas pressure and cluster mass. As a result, the sample of sources 
with estimates of the ambient gas pressure and cluster mass is rather small. 

Results for radio-loud quasars and radio galaxies are analyzed separately when there are 
enough data points, as is the case in the study of the beam power. As noted above, the number of 
sources with thermal pressure and cluster mass estimates is rather small (16 lobes with detections 
and 14 with bounds). Thus it is not practical to separate quasars and galaxies in the study of 
these parameters. 

Cygnus A (3C 405) is the only source with a redshift ~ in our sample, and in some cases, 
results seem to depend on whether or not this source is included. Thus, subsamples both with and 
without Cygnus A are analyzed. 

3. Luminosity in Directed Kinetic Energy of the Jet 

3.1. Theory 

The jet or beam of an FRII radio source refers to the collimated outflow which carries energy 
from the AGN. When the jet impinges upon the external medium, a strong shock is formed and 
the kinetic energy of the jet is deposited in the vicinity of the shock front, which is marked by 
the radio hot spot. In principle, one can calculate the beam power, dE/dt, carried by the jet 
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by studying the propagation of the hot spots. However, in practice it is difficult to use hot spot 
properties to estimate the beam power, since the hot spot is generally not resolved, and the 
position and properties of radio hot spots vary on short time scales (Laing 1989; Carilli, Perley, 
and Dreher 1988; Black et al. 1992.). 

The model given by (Daly 1994, 1995) bypasses these difficulties by using properties of the 
more stable radio lobes to estimated the luminosity in directed kinetic energy, or the beam power, 
Lj. Observations and numerical simulations indicate that radio lobes of powerful extended radio 
sources propagate supersonically relative to the ambient medium (Alexander & Leahy 1987; 
Prestage & Peacock 1988; Cox, Gull, & Scheuer 1991; Daly 1994; WDW97a,b). Leahy (1990) 
shows that the properties of the radio lobe can be used to estimate the rate of energy input. His 
equation may be rewritten as: 

Lj = ^a\v L P L . (1) 

Here cll is the cross-sectional radius of the radio lobe, vl is the lobe propagation velocity, and the 
FRII source is assumed to be cylindrically symmetric about the radio axis. Using typical units, 
Lj can be expressed as 

Lj = 3.6 x 10 44 erg ^ l ( ^-Y (-) ( Tn~r~~ 2] > (2) 

3 6 Vkpc/ V c J Vl0- 10 dyne cm" 2 / w 

All three parameters used to estimate Lj in eq. (0) can be estimated from radio data (see, 
for example, Daly 1995, or WDW97a,b). The lobe radius ai can be readily measured from 
the radio map (see WDW97a,b for detail). The lobe propagation velocity vl can be estimated 
using multi-frequency observations of radio lobes and bridges, by modeling the spectral aging of 
relativistic electrons due to synchrotron and inverse Compton losses (c/. WDW97a,b; Wan & Daly 
1998a,b; Myers & Spangler 1985; Alexander & Leahy 1987; LMS89; LPR92; Carilli et al. 1991). 
The lobe pressure Pl is given by 

P L = Q r i. 5 + ^^ L) (3) 

where B m i n is the magnetic field strength in the radio lobe under the minimum-energy condition 
(c/. Miley 1980; Pacholczyk 1970), b is the ratio of the true magnetic field strength B to the 
minimum-energy magnetic field strength: B = bB m i n , and the equation is in cgs units. WDW97b 
estimate b to be about 0.25, with a source-to-source dispersion less than about 15%. This value 
of b is consistent with values obtained by other authors (c/. Carilli et al. 1991; Perley & Taylor 
1991; Feigelson et al. 1995), and is used throughout this study. 

The radio spectral index of an FRII source is used in the spectral aging analysis to estimated 
the lobe propagation velocity vl- It is shown in WDW97b that there is a clear correlation between 
the radio spectral indices of the sources in this sample and redshift. The spectral index a can 
be expressed as a oc (1 + z) s . The best fitted value of s is 0.8 ± 0.2 with a reduced \ 2 of 1.7. 
Such a correlation could be caused by intrinsic curvature in the initial electron energy spectrum 
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or effects of inverse Compton cooling on the hot spot spectral index (see WDW97b). Thus, the 
observed spectral index may not be the appropriate index to used in spectral aging analysis. This 
lead WDW97b to also consider correcting the observed spectral indices to zero redshift using the 
empirical correlation between a and (1 + z) for the sample. 

The beam power Lj can be used to estimate the total time that the AGN will be active, and 
will produce highly collimated outflows. Following Daly (1994, 1995) and Guerra & Daly (1996, 
1998), the total lifetime of an outflow, defined as t*, is related to the energy extraction rate Lj by 
a power law: 

i* = CLjM 3 , (4) 

where /3* is a parameter that has implications for models of energy extraction from the central 
engine. The current best estimate of /?* is 1.75 ±0.25 (Guerra, Daly, & Wan 1998). (Note that the 

(3* used here is not related to (3q used in §|| and §[| for the King density profile.) The characteristic 

—8 /3 

size of an FRII source is related to its lifetime as D± = v^t* = CvlL- * , where C is the 
normalization in eq. (0). This normalization factor is chosen so that at z ~ 0, the characteristic 
size of Cygnus A matches the observed average lobe-lobe size for sources at this redshift. That is, 
(v£)CL- for one side of Cygnus A is added to (vl)CL- on the other side of Cygnus A, 
and this is set equal to 2(< D > 2 =o) = 2 x (68 ± 14) kpc. This equation is then solved for the 
normalization factor C; the uncertainty on C is indicated in Table 1. This normalization is then 
used in equation (Q) to obtain an estimate of the total time that the AGN will be producing 
large-scale jets. The value for t* estimated using independent information from each side of the 
source should of course be equal, and are generally very close in value. The value adopted for the 
beam power, Lj, of Cygnus A used above is listed in Table 1, x>l for Cygnus A is obtained from 
Table 1 from WDW97b, and a value of < D > z=0 = (68 ± 14) kpc is adopted from Table 3 of 
Guerra, Daly, & Wan (1998). 



3.2. Empirical Results 

The beam power, Lj, has been estimated using eq. @ for the lobe propagation velocity 
estimated using the observed radio spectral index, and the redshift-corrected radio spectral index; 
are both listed in Table |]. Typical values of Lj are 10 45 h~ 2 erg s . 

No strong correlation is found between the beam power Lj and the linear size of the source, 
represented by the core-hot spot separation r. Figure || is a log-log plot of Lj vs. r, where the 
observed radio spectral index a is used. The figure of Lj vs. r after applying the redshift-correction 
on a is almost identical to Figure ^, and is thus not shown here. 

For radio-loud quasars, the best-fitted line in Fig. || has a slope of about zero. A Spearman 
rank correlation analysis shows that the correlation coefficient between Lj and r is —0.03(7.72%) 
when no redshift-correction on the spectral index is applied, and —0.01(1.78%) when the spectral 
index correction is applied. Here the number in parenthesis is the significance of the correlation. 
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These results suggest that the correlation between Lj and r for radio-loud quasars is insignificant, 
either with or without the a — z correction. That is, Lj is independent of r for the radio-loud 
quasars. 

For the radio galaxies, Lj appears to increase slightly with r in Fig. ||, and results from the 
Spearman rank correlation analysis hint a marginally significant correlation between Lj and r. 
The correlation coefficient is about 0.3, with a significance level of about 90%, which holds with 
or without the a — z correction. However, note that the best-fitted line for radio galaxies in Fig. || 
has a slope only 2a away from zero, and that a correlation coefficient of about 0.3 is rather weak. 
Thus it appears that Lj is at most weakly dependent on r for the radio galaxies. These results 
for radio-loud quasars and radio galaxies are consistent with the assumption that Lj is roughly 
constant over a source's lifetime. 

The relationship between Lj and redshift and radio power (P r ) has been studied in detail 
in Wan &; Daly (1998a). Increases of Lj with z and with P r are observed (see Figures || and 
^). However, the radio power of the sources in the sample also increases with redshift since the 
sources used in this study all come from the flux-limited 3CR survey. As a result, which of the 
Lj — z and Lj — P r correlations is more significant, and the role of radio power selection effects, 
needs to be considered carefully. Wan & Daly (1998a) used two parameter fitting and partial 
rank analysis to investigate this. The two parameter fit Lj cx Pyjg(l + z) nz yields the following n z 
values: n z = 1.45 ± 0.32(0.60) for all galaxies, n z = 3.83 ± 0.93(1.73) for all galaxies except for 
Cygnus A, and n z = 9.56 ± 1.37(3.91) for all radio-loud quasars. 

The number in parenthesis is the uncertainty on n z times \/rx 2 , which includes the effect of 
a reduced x 2 ( r X 2 ) that is greater than one for the fit. This notation will be used throughout this 
paper. 

The values of n z are more that 2a away from zero, for both radio galaxies and radio-loud 
quasars. This means that the Lj — z correlation is unlikely to be caused purely by radio power 
selection effects. That is, a real increase of Lj with redshift exists, though the magnitude of this 
redshift evolution is not well determined. 

The Lj — P r correlation, on the other hand, is rather poorly constrained. Results from the 
one parameter fit suggest Lj goes roughly proportional to Pi 73 (Figure ^). The two parameter fit 
Lj oc Pi7g(l + z) nz , which takes into account the redshift dependence of radio power, gives values 
of n p with large uncertainties: n p = 0.76 ± 0.11(0.20) for all galaxies, n p = 0.25 ± 0.22(0.41) for all 
galaxies except for Cygnus A, and n p = —0.54 ± 0.29(0.82) for all radio-loud quasars. Thus, the 
exact relationship between Lj and radio power is not clear at the present time. Note, however, 
large amounts of scatter are seen in the Lj vs. Pi 78 plot (Figure |^), which is also reflected in the 
large reduced x 2 °f t ne one parameter fit. For the same beam power, the radio power can vary 
by about a factor of ten. This suggests that radio power is not an accurate measure of the beam 
power. 

The fact that Lj increases with redshift implies that the total lifetime of an outflow, t+, 
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estimated using eq. (||), decreases with redshift. The values of the full lifetime of the source 
are listed in Table [lj. Typical values of t* are about 10 7 years. Figure |6| plots t* as a function 
of redshift. It is clear that decreases with z, which can explain the fact that the average size 
of powerful extended 3CR sources decreases with redshift for z > 0.3 (see Guerra & Daly 1998), 
while the lobe propagation velocity vl increases with redshift. Clearly if vl increases with redshift, 
and the mean source size decreases with redshift, then t* must decrease with redshift since 
D* ~ VLt*. The sources at high redshift produces more powerful jets (with larger Lj) for a shorter 
period of time (with smaller i*), which results in smaller average sizes than low-redshift sources. 



4. Thermal Pressure of the Ambient Gas 
4.1. Theory 

The studies of WDW97a and b show that both the ambient gas density and temperature of 
an FRII source can be estimated from radio data. The thermal pressure of the ambient gas Pth is 
obviously proportional to the product of the density n a and the temperature T. Interestingly, as 
shown below, although the lobe propagation velocity enters into n a and T, it cancels out in the 
product, so the thermal pressure of the ambient gas can be estimated using single frequency radio 
data if the radio emission from the bridge is mapped over a large enough region of the bridge. This 
offers a completely new method to estimate the thermal pressure of the ambient gas surrounding 
high-redshift radio sources using only single frequency radio data. 

It can be complicated and time consuming to estimate the ambient gas density in the vicinity 
of classical double radio sources using X-ray observations, especially for sources at high redshift. 
An attractive alternative method of estimating n a in the vicinity of a radio source is to use the 
radio data. Given that the radio lobe represents a strong shock front, the lobe pressure Pl, the 
ambient gas density p a , and the lobe propagation velocity are related: the strong shock jump 
conditions imply that Pl ~ 0.75p a v L holds, where Pl is the non-thermal pressure inside the radio 
lobe, and p a is the mass density of the ambient gas. The electron number density of the ambient 
gas n a can be estimated using 

Pl (5) 
\Apm v v\^ 

where p is the mean molecular weight of the gas in AMU, with a value of 0.63 when solar 
abundances are assumed, and m p is the proton mass. 

The ambient gas temperature can be estimated by careful studies of the radio bridge, as 
discussed in detail in WDW97a. A brief description of the theory is as follows. When the 
non-thermal pressure of the radio bridge is much greater than the ambient gas pressure, the bridge 
undergoes supersonic expansion. This lateral expansion can be treated as a blast wave and the 
expansion velocity is determined by ram pressure confinement. This causes the width of the bridge 
to vary with time approximately as i 1 / 2 (e.g. Begelman & Cioffi 1989; Daly 1990; WDW97a). 
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Since the lobe propagation velocity vl is roughly constant during a given source's lifetime {e.g. 
LPR92; WDW97a,b), the lobe front will advance a distance of vjj, during t as the bridge expands. 
This means that at a distance x from the hot spot, the bridge expansion time is t = x/vl- Thus 
the shape of the radio bridge roughly follows a square root law, with the bridge width a(x) oc x 1 ! 2 . 
This behavior was in fact found by WDW97a. As the bridge expands, the pressure inside decreases 
until it becomes comparable to the ambient gas pressure. When this occurs, the lateral expansion 
velocity becomes sonic, causing a break in the functional form of a{x). The width of the bridge 
starts to deviate from the square root law and becomes roughly constant as equilibrium is being 
reached. At the point where such a break occurs, the lateral expansion velocity is approximately 
sonic. Thus the sound speed of the ambient gas (c s ) can be estimated using 
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where v* 



■j s \ h denotes the lateral expansion velocity at the break. 

The Mach number of lobe propagation is defined as M = vl/c s . Given the expression for c 
in eq. @, the Mach number of lobe propagation is: 

da x - 1 
dx 



M 



VL_ 
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a b '' 



(7) 



where a b is the width of the bridge at the break, and x b is the position of the break relative to the 
hot spot (see WDW97a). It can be seen from eq. (j^) that M depends only on the geometrical 
shape of the radio bridge. 



The ambient gas temperature T can be expressed as: 



fim p 2 
-c. cx 



da 
dx 



v\ cx v 2 L M-\ 



(8) 



where 7 is the ratio of specific heats of the ambient gas. 

Multi- frequency radio maps are required in order to use eqs. and (||) to estimate the 
ambient gas density and temperature of an FRII source, since both n a and T depend on wj,. 
However, there is no such requirement in order to estimate the thermal pressure of the ambient 
gas. The ambient gas pressure in the vicinity of the radio lobe is simply 



th 



(n a + rii)kT oc ( 



v 2 T M- 2 o, Pl 



M 2 ' 



(9) 



with Pl estimated using eq. (H) and M estimated using eq. (0). Here rii is the number densities of 
ions, and for a gas with solar abundances, the electron density is n a = 1.21m. Substituting in the 
numerical constants in eqs. (||) and (||), the thermal pressure is given by 



th 



0.8 



Pl_ 

M 2 ' 



(10) 
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for a gas with a specific heat ratio 7 = 5/3. 

The thermal pressure estimated using eq. ([R]) depends only on the lobe pressure Pl and the 
geometrically determined Mach number M, both of which can be obtained from single frequency 
observations of the radio lobe and bridge. No spectral aging analysis is needed in order to estimate 
Pth, since the dependences of n a and T on vl cancel. As a result, P t h is also independent of 
whether or not the redshift-correction on the radio spectral index is applied. 

At low-redshift, X-ray data can often be used to estimate n a and T, and hence Pth- However, 
high-resolution X-ray data are often difficult to obtain for sources at high redshift. This new 
method of estimating Pth, using single frequency radio data, offers an attractive alternative. It 
provides a powerful tool to study the environments of powerful classical double radio sources, 
especially those at high-redshift. Since it has been shown that these sources are in cluster-like 
gaseous environments, the sources may be used to study evolution of gas in clusters of galaxies (as 
discussed by Daly 1995; WDW97a,b). 

4.2. Empirical Results 

The value of Pth, the thermal pressure of the ambient gas in the vicinity of the radio 
lobe, is listed in Table [fl. Most sources in the sample have Pth on the order of (10 -11 to 
10 -10 ) /1 4 / 7 dyne cm -2 , which is typical of gas in low-redshift clusters of galaxies. This is consistent 
with results obtained by WDW97(a,b), who find cluster-like density and temperature for the 
ambient gas of the FRII sources in the sample. Note that the thermal pressures obtained here do 
not depend on a spectral aging analysis, whereas density and temperature estimates do. 

The X-cluster around Cygnus has been observed by numerous authors (Arnaud et. al. 1984; 
Carilli, Perley, & Harris 1994). The thermal pressure of the ambient gas near its radio lobe is 
estimated to be about 10~ 10 dyne cm -2 (see Carilli, Perley, & Harris 1994). This estimate from 
X-ray data is consistent with the thermal pressure estimated here for Cygnus A using the new 
method. 

Figure || plots P t h as a function of the core-hot spot separation r. It can be seen that P t h 
decreases with r. A negative pressure gradient is expected for an isothermal gas distribution that 
follows the King density profile, which is suggested to be the ambient gas distribution for sources 
in our sample (WDW97a,b). For this gas distribution, the thermal pressure decreases with r as 

-§A> 

(11) 

where Pth, c is the thermal pressure at the center of the cluster, and r c is cluster core radius. 

The pressure gradient seen Figure || appears to be consistent with that expected for an 
cluster gas distribution that is isothermal and follows a King density profile, with a cluster core 



Pth(r) = Pt 



th, c 
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radius r c of ~ (50 to 150) h 1 kpc. The sources in our sample have values of r ranging from 
(25 to 250) h~ l kpc. Within this radius range, the pressure profile given by eq. ( |TT| ) has an 
average slope of (~— 1.4 to — 0.7), for r c values of (50 to 150 h~ x kpc), assuming /3 = 2/3. These 
expected slopes are consistent with the best-fit slopes of ~ — 1.0 ± 0.3 in Figure []. 

There is some hint of a redshift evolution of the cluster core radius from the data. Figure [6| 
plots Pth vs. r in two redshift bins, z < 0.9 and z > 0.9. The division at z = 0.9 is chosen so that 
the two redshift bins cover about the same range in redshift and contain about the same number 
of data points. It appears that Pth decreases less rapidly with r at high redshift than at low 
redshift, suggesting a larger r c at high redshift. The best-fit slope of ~— 1.4 in the low-redshift bin 
is consistent with r c ~ 50/i _1 kpc for the pressure profile given by eq. (pd|), whereas the best-fit 
slope of about —0.4 in the high-redshift bin is consistent with r c ~ 250 h -1 kpc. This result is still 
preliminary since only sources with detections of Pth are included in the fits. Better estimates of 
Pth for sources currently with only upper bounds on P t h will help to better determine whether r c 
is evolving with redshift. 

The results obtained above are consistent with results obtained by WDW97b. They study the 
ambient gas density profile and find that r c ~ 50 /i" 1 kpc when high- and low-redshift source are 
considered together. They also find that that data are consistent with a constant core gas mass 
model where r c increases with redshift roughly as r c oc (1 + z) . Note that results from the study 
of Pth do not depend on spectral aging analysis, whereas those from the study of n a do. It is thus 
encouraging that consistent results are obtained from the two studies. 

For the following analysis, we use both a non-evolving core radius of 50 /i" 1 kpc, and an 
evolving core radius of 50 (1 + z) 1 ' 6 /! -1 kpc, with a focus on the latter. 

The thermal pressure at the center of the cluster (Pth, c ) can be estimated by scaling the 
pressure in the vicinity of the radio lobe (Pth) to the cluster center. This follows because the 
studies of Daly (1995) and WDW97a,b indicate that these radio sources are located at the centers 
of clusters of galaxies, so the core-hot spot separation can be used as an estimate of distance from 
the cluster center. The current data on T, n a , and Pth are all consistent with the gas distribution 
being isothermal and following a King profile. Thus we use the pressure profile for such a gas 



distribution, as given by eq. (11), to estimate the central thermal pressure. The central pressure is 
simply 



th, c 



Pth(r) 



1+ ( - 



r .2" 



(12) 



The values of Pth.c for the sources in our sample are listed in Table |], where a value of (5q = 2/3 
has been assumed. Most sources have values of Pth,c around 10 -10 /i 4 / 7 dyne cm" 2 , which is rather 
typical of gas pressure in the core regions of low-redshift galaxies clusters. 

Results on the redshift evolution of the central gas pressure Pth,c are somewhat uncertain 
because Pth,c seems to correlate with both redshift and radio power when either correlation is 
considered separately (see Figures and 0) , and it is not clear at present which correlation is more 
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significant. For the non-evolving core radius model, a two-parameter fit of Pth c oc (1 + z) nz Pyfg 
yields n z = 0.29 ± 0.49(0.67) and n p = 0.66 ± 0.19(0.26) when all sources are included, and 
n z = 1.51 ± 2.11(2.99) and n p = 0.34 ± 0.58(0.82) when Cygnus A is not included. For the evolving 
core radius model, the two-parameter fit gives n z = —1.22 ± 0.49(0.65) and n p = 0.88 ± 0.19(0.25) 
when all sources are included, and n z = 1.56 ± 2.11(2.82) and n v = 0.15 ± 0.58(0.78) when 
Cygnus A is not included. Note that Cygnus A appears to lie in a cooling-flow region while most 
of the other sources are not in cooling- flow regions (see WDW97a; §|5|) . Thus fits with and without 
Cygnus A are both performed. In any case, the large uncertainties on n z and n p make it hard to 
determine the magnitude of the redshift evolution of Pth,c- 

The thermal pressure can be used to predict the amount of cosmic microwave background 
(CMB) diminution, also known as the Sunyaev-Zel'dovich (S-Z) effect, that is expected from the 
cluster. The reduction in the cosmic radiation (AI U ) in the direction of a cluster is given by (e.g. 
Sunyaev & Zel'dovich 1980; Sarazin 1988; Rephaeli 1995) 

AI U /I U = G{u)y. (13) 

Here I u is the specific intensity of the CMB at the observing frequency v, y is the Comptonization 
parameter, and the function G(oj) is 



G{u)) 



it! 



e^ + l 
e u - 1 



(14) 



- 1 

where u = hv/kT r , and T r = 2.73 K is the CMB temperature (Mather et al. 1990). 
The Compton y parameter is given by 

y = [ -^ra T n e dl = — ^ ( — — — ) I ' Pthdl. (15) 

J m P cr m P c z \n e + n 



where ot = (87r/3)[e 2 /(m e c 2 )] 2 is the Thompson electron scattering cross section, T is the cluster 
gas temperature, n e is the electron density, and is the ion density. Using the thermal pressure 
in a cluster with an isothermal-King gas distribution (eq. [|ll|), this gives 



th, c 



r(3/? /2 - 1/2) 



y « 6.3 x 10- 5 ! (1 + s 2 ) -3A,/2+ 1 /2 j (16a) 

\10- 10 dyne cm" 2 / \0.25Mpc/ r(3/3 /2) 

where x = r/r c , and the gas is taken to have solar abundance so that n e = 1.21nj. For a typical 
value of = 2/3, 

y « 1.12 x 10- 4 ( — ^ ( (1 + x 2 )' 1 / 2 . (16b) 

y VlO" 10 dyne cm" 2 / V 0.25 Mpc/ v ' y ' 



At low frequency, clusters with strong radio sources are generally avoided for measurements 
of the Sunyaev-Zel'dovich effect because emission from the radio source masks the microwave 
diminution. Such contaminations from radio sources are reduced at high frequency since emission 
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from steep-spectrum radio sources, such as the FRII sources studied here, decreases rapidly with 
increasing frequency. The Sunyaev-Zel'dovich Infrared Experiment (SuZIE) can measure the S-Z 
effect around 140 GHz, where the amount of CMB intensity diminution is near its peak. At this 
frequency, a cluster with an isothermal-King density profile will cause a CMB intensity diminution 
of about 



Pth, 



AJ„ » 7.0 x 10~ 19 erg s^Hz^OTrtr -1 ( ^ - , . 

& V 10~ 10 dyne cm" 2 J \ 0.25 Mpc 

v r (3/3o/2 - 1/2) ^ + x2) _ 3/3()/2+1/ ^ 



(17) 



r(3A>/2) 

and for a typical value of [3q = 2/3, 

A/, « 1.24 x 10- 18 erg s^Hz-W 2 **- 1 ( T7r ^ ^ f ^— ) (1 + x 2 )' 1 ' 2 . (18) 

Vl0- 10 dyne cm" 2 / V0.25Mp C y v ; v ; 

For a detector with a FWHM beam size 9b, the total amount of CMB diminution within the 
beam, defined as AF U , is roughly 

AF„ « / • 24mJy f ^ V ( S) ( — , (19) 

J y U-7 arcmin/ VlO" 10 dyne cnr 2 / V 0.25 Mpc/ ' V ; 

where / is the beam dilution factor, defined as the ratio of the average AI U within the beam to 
the peak value at the cluster center; and the SuZIE FWHM beam size at 140 GHz of 1.7 arcmin 
(e.g. Holzapfel et al. 1997) is used to calculate the numerical value. For a Gaussian beam, when 
the HWHM beam size corresponds to lr c , / ~ 0.7, and when the HWHM beam size corresponds 
to 2r c , / 0.5 (see Rephaeli 1987). At z < 2, the SuZIE beam radius is less than ~ 320/i _1 kpc 
for go = with no cosmological constant. Thus a cluster with Pth,c ~ 10~ 10 dyne cm -2 and 
r c ~ 0.25 Mpc will have a AF U of about (15 to 20) mJy within the SuZIE beam. The clusters 
surrounding the FRII sources in our sample have Pth,c ° n the order of 10 -10 /i 4 / 7 dyne cm -2 and 
the data also suggest that r c increases with redshift, reaching about 300 h~ l kpc at z~2. The 
expected CMB diminution within the SuZIE beam at 140 GHz for these clusters ranges from 
several to tens of mJy. These clusters make good candidates for SuZIE observations if the fluxes 
from the radio sources at 140 GHz or their uncertainties are small compared with the expected 
S-Z effect signals. 

We are currently in the process of searching for high frequency data on the radio sources 
in our sample in oder to identify possible SuZIE observation candidates. One likely candidate 
that comes from a preliminary search in the published literature is 3C239. The expected CMB 
diminution for the cluster surrounding it is about 30 h~ 3 ^ 7 mJy within the SuZIE beam at 140 
GHz. In contrast, extrapolation of the observed spectrum of 3C239 as given by observations at 
178 MHz, 10.7 GHz, and 14.9 GHz (Kellermann k Pauliny-Toth 1973; Genzel et al. 1976; Laing, 
Riley, and Longair 1983), gives a 140 GHz flux of only 6 mJy. This is likely to be an overestimate 
since synchrotron aging and inverse Compton cooling can both cause the high frequency spectral 
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index to become steeper than that at low frequency. The effects of inverse Compton cooling can 
be especially important since this source is at high redshift (z = 1.79), where the energy density 
of the microwave background is much higher than at low redshift. Thus it appears that emission 
from this radio source is weak compared with the expected CMB diminution. Observations of the 
radio source at more frequencies above 14.9 GHz can help to better constrain its 140 GHz flux. 



5. Gravitational Mass of the Host Cluster 
5.1. Theory 

The powerful classical double radio sources studied here are in cluster-like gaseous 
environments (see WDW97a, b). The studies of these sources provide information on density 
and temperature of the ambient gas. Thus it is possible to estimate one of the most important 
parameters of the cluster, the total gravitational mass, including dark matter. 

The total mass can be estimated using the density and temperature profile of the intracluster 
medium (ICM) if the gas is in hydrostatic equilibrium. The sound crossing time in the ICM is 
usually short compared to the age of a high-redshift cluster, and the morphology of the X-ray 
emission from many low-redshift clusters is often smooth. Thus it is generally believed that 
hydrostatic equilibrium is a good approximation of the state of the ICM for many clusters that 
are not cooling flow clusters. Assuming spherical symmetry, hydrostatic equilibrium requires 

1 dP GMAr) . . 

-— = 20 

p g dr r z 



where p g is the gas density and Mj(r) is the total gravitating mass within r. This means that 

(21) 



», ^ kT ( r ) 
M t (r) - 



Gprrip 



d\ogp g {r) d log T(r) 



d log r d log r 



where k is the Boltzman constant, p is the mean molecular weight of the gas in amu, m p is the 
proton mass, and T(r) and p g (r) are the temperature and density profiles of the cluster gas. 

The most commonly used hydrostatic model of the ICM is the isothermal /3-model (Cavaliere 
& Fusco-Femiano, 1976,1978; Sarazin & Bahcall 1977), where the clusters gas is isothermal and 
the density profile follows a modified King model, i.e., n a = n c [l + (r/r c ) 2 ] -3//2 ^°. Here n c is the 
core density and r c is core radius of the gas distribution. Using this model for the ICM, eq. (gjj) 
becomes 

JH W = J**.. r . r .^^. (22) 
&pm p 1 + (r/r c y 

Results from numerical simulations of cluster formation suggest that departures of the ICM from 
hydrostatic equilibrium are usually small, and the mass estimated using the standard /3-model is 
rather accurate {e.g. Navarro, Frenk & White 1995; Schindler 1996; Evrard, Metzler, & Navarro 
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1996). Note however, significant temperature decline is observed to occur at outer regions of 
clusters (e.g. Markevitch et al. 1997). Thus, we only use eq. (^) to estimate the cluster mass 
within r rather than extrapolating to large radii. 



The mass estimated using eq. ( 22 ) is not accurate if the cluster is cooling at the point where 
the cluster temperature is measured. Inside a cooling flow region, the temperature measured 
does not indicate the gravitational potential, since hydrostatic equilibrium conditions do not 
apply. Considering cooling only by thermal bremsstrahlung, which dominates other mechanisms 
at typical clusters temperature, the cooling time of the cluster gas can be estimated using 



t cool « 2.1 x 1(T 2 ( ) ( ) Gyr, (23) 



( T V /2 ( »« r 1 

U0 7 Ky Urn" 3 , 

where T and n a are the ambient gas temperature and electron density, respectively (see WDW97a). 
Note this equation does not depend on the Hubble's constant, provided that the estimates of n a 
and T are not dependent on the Hubble's constant. 

Given the ambient gas temperature and density estimates, the cooling time for the sources in 
our sample can be estimated (see WDW97a). A cooling time less than the age of the universe at 
that redshift means that the FRII source is in a cooling flow region. To estimate the age of the 
universe at a given redshift, an open empty universe is assumed, with a current age of 14 Gyr, 
which corresponds to h ~ 0.7 for qo = with zero cosmological constant. 

Most sources in our sample do not appear to be cooling at the position where the temperature 
is measured. Thus their mass estimates are likely to be reliable. A few sources, including the 
low-redshift source Cygnus A, appear to lie in cooling flow regions. 

Figures || plots the total cluster mass, including dark matter, within radius r, Mt(r), as a 
function of the core-hot spot separation r. Since the sources appear to lie at the centers of clusters 
of galaxies (Daly 1995; WDW97a,b), the core-hot spot separation is used to estimate the distance 
from the cluster center. The radio spectral index is redshift-corrected and the cluster core radius 
is taken to be r c = 50(1 + z) 16 h^ 1 kpc (see §4.2). Clusters that are cooling at the position where 



the temperature is measured are marked with pentagons. The plots of Mt(r) vs. r for other 
models, either with or without an a — z correction and/or a core radius evolution, are presented in 
Figures ^ through ||. 

The cluster mass increases with r for all models. The increase is more rapid for models with 
an evolving core radius than models with a fixed core radius. This merely reflects the fact that 
with an increasing r c , more sources lie within the core region, where the increase of mass with 
radius is rapid, and varies roughly as Mar 3 . 

To determine the redshift evolution of the cluster mass, a three parameter fit of 
Mf(r) oc r ni (l + z) n2 P^j\ is performed since the cluster mass within r, Mt(r), is a function of r, 
and may also be affected by radio power selection effects. Results are listed in Table [2] for all 
models considered. A two parameter fit of Mj(r) oc r ni (l + z) 112 is also performed, and results are 
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listed in Table ||. For completeness, fits including all sources, all sources except Cygnus A, and 
all sources that are not in cooling flow regions, are all listed in the table. Note that Cygnus A 
appears to be in a cooling flow region for all the models considered. Results obtained excluding 
cooling flow clusters are probably the correct ones to consider. 

The sample of clusters with mass estimates is rather small, especially when only non-cooling 
flow clusters are considered. Thus it is not surprising to see that the best-fit values of n<i have 
rather large uncertainties, which makes it hard to draw definitive conclusions about the redshift 
evolution of the cluster mass. Note though, the current data do not indicate any negative evolution 
of the cluster mass out to a redshift of about two, contrary to what is expected in cosmological 
models with a high density parameter Q m . This is consistent with results obtained by this group 
(Daly 1994; Guerra & Daly 1996, 1998; Guerra, Daly, & Wan 1998; and Daly, Guerra, & Wan 
1998), who study the characteristic size of powerful classical double radio sources and find the data 
suggest a low value of Q m , which is greater than 3a (check) away from a universe with Q m = 1 
(see Guerra, Daly, & Wan 1998). 

There are some indications from the data that the redshift evolution of the cluster core radius, 
and the redshift-correction on the radio spectral index are favored. The gas mass within r, defined 
as M g (r), can be estimated using 

M g (r) = Airp c rl [x - tan -1 (x)] , (24) 

where p c is the central gas density, r c is the core radius, x = r/r c , and a typical value of 0o = 2/3 
is used for the King density profile. Knowing the gas mass and the total mass within r, the 
gas mass fraction within r can be estimated. Figure || plots the gas mass fraction within r vs. 
redshift, where the radio spectral index is redshift-corrected and the core radius is taken to be 
r c = 50(1 + z) ls hT 1 kpc. 

The values of gas fraction shown in the figure are consistent with observed values for 
inner regions of many clusters (e.g. Donahue 1996), whereas those for models without a 
redshift-correction on the radio spectral index, and without a core radius evolution seem to be low 
compared with observed values. 

Several factors may contribute to the slow decrease of gas fraction with redshift that is seen 
in Figure ||. First, the increase of cluster core radius with redshift means that the gas fraction 
estimated for the high-redshift clusters is over a larger fraction of the cluster core than for the 
low-redshift clusters. The gas fraction estimated at high redshift mainly represents the gas fraction 
inside the cluster core, and it is known that gas fraction increases with increasing radius in clusters 
(e.g. David, Jones, & Forman 1995). Second, the mass contribution from galaxies is higher at 
the cluster center than at large radii (e.g. Loewenstein & Mushotzky 1996). Thus by sampling a 
larger fraction of the cluster core at higher redshift, a larger fraction of the total baryon mass is 
not taken into account at higher redshift. Further, the decrease in gas fraction with redshift could 
be due to the stripping of gas from galaxies in the cluster core over time. Futhermore, the cluster 
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gas becomes more concentrated toward the cluster center as it cools, which also causes the gas 
fraction in the cluster center to increase with time. 

These results on cluster mass and gas fraction obtained above are still preliminary due to the 
small size of the sample. More sources with estimates on the ambient gas temperature, and hence 
cluster mass, will help to test these results. 

6. Summary and Discussion 

Several key parameters of an FRII source and its gaseous environment are studied in this 
paper. 

Direct estimates of the beam power, which measures the energy extraction rate from the AGN 
by the jet, are obtained. Typical beam powers of about 10 45 h~ 2 erg s _1 are found for the sources 
in the sample. No strong correlation is seen between the beam power and the linear size of the 
source, which is consistent with the beam power being roughly constant throughout the lifetime 
of a source. 

There is a trend for the beam power to increase with redshift, which is significant even after 
excluding radio power selection effects. The magnitude of this redshift evolution, however, is 
not well determined by the current data. It is well known that the quasar luminosity function 
undergoes strong evolution between z ~ and z ~ 2, with the high-redshift quasars being more 
luminous than their low-redshift counterparts (e.g. Schmidt & Green 1983; Yee & Ellingson 1993; 
La Franca & Cristiani 1997). This suggests that high-redshift AGNs are more powerful than low 
redshift ones. Thus it is perhaps not surprising to find that the beam power, or energy per unit 
time channeled into the jet by the AGN, is also higher at high redshift. 

The relationship between the beam power and the radio power is not well constrained after 
their correlations with redshift are taken into account. The two parameter fit of Lj oc (1 + z) Uz P^f 8 
yields values of n v with large uncertainties. Thus it is not clear at present how the beam power and 
the radio power are related. However, the large amount of scatter seen in the beam power-radio 
power relation suggests that radio power is not an accurate measure of beam power. 

The beam power Lj can be used to estimate the total lifetime of an outflow produced by 

an AGN. Following Daly (1994) and Guerra & Daly (1996, 1998), the total time for which the 

—3 /3 

outflow will occur, t*, is related to Lf t+ oc L- , where /3* is estimated to be about 1.75 ± 0.25 
(Guerra, Daly, & Wan 1998). Typical lifetimes of about 10 7 years are obtained using this relation 
for the sources studied here. The lifetime of the outflow decreases with redshift, which explains 
the decrease of the average size of powerful extended 3CR sources with redshift. The sources at 
high redshift produces jets for a shorter period of time, which results in smaller average sizes than 
low-redshift sources. 

A new method of estimating the thermal pressure of the ambient gas in the vicinity of powerful 
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classical double radio source using only single frequency radio data is presented. The pressure is 
given by the product of the ambient gas density, estimated using ram pressure confinement of the 
radio lobe, and the ambient gas temperature, estimated using the Mach number for the source 
and the lobe advance speed. It turns out that the lobe propagation velocity cancels out of the 
product n a T, and the ambient gas pressure can be estimated by studying the shape of the radio 
bridge, and the non-thermal pressure of the radio lobe (see Thus, the thermal pressure, the 
product of the ambient gas density and temperature, depends only on the non-thermal pressure in 
the radio lobe and the geometrically determined Mach number of lobe advance, both of which can 
be estimated using single frequency radio data. This new method to estimate the thermal pressure 
does not require a spectral aging analysis, and provides a powerful tool to probe the environments 
of FRII sources. The thermal pressure estimated for Cygnus A using this new method agrees with 
that obtained using X-ray data. 

Thermal pressures on the order of 1CP 10 /i 4 / 7 dyne cm -2 , typical of gas in low-redshift clusters 
of galaxies, are found for the gaseous environments of the FRII sources studied here. There are 
hints from the current data that the gradient of the composite pressure profile is less steep at high 
redshift than at low redshift, which can be explained by an increase of the cluster core radius with 
redshift. The current data are consistent with a core radius evolution from ~50/i _1 kpc at z~0 to 
~250/i~ 1 kpc at z~2, which agrees with results obtained from studies of the ambient gas density 
(WDW97b). WDW97b find that the data can be described by a model where the core gas density 
decreases and the core radius increases so that the core gas mass remains roughly constant. 

The thermal pressures obtained here can be used to estimate the the amount of 
CMB diminution expected from the clusters surrounding the FRII sources in the sample. 
Contaminations from the radio sources can be reduced by observing at high frequency, such as 
140 GHz, since emission from the radio sources decreases rapidly with increasing frequency. The 
cluster surrounding a source in our sample can be detected by SuZIE observations at 140 GHz 
if the flux from the radio source at this frequency or its uncertainty is small compared with the 
expected S-Z effect signal. A search for high-frequency data of the sources in our sample is ongoing 
in oder to identify possible SuZIE observation candidates. Preliminary results suggest 3C239 to 
be a likely candidate. 

The gravitational or total mass of the surrounding cluster is estimated for the sources in the 
sample, assuming hydrostatic equilibrium conditions for the gas. Masses of up to 10 14 M Q are 
found for the central regions (r 250 /i -1 kpc) of the clusters, consistent with typical values for 
low-redshift clusters. The redshift evolution of the cluster mass is not well determined. Current 
data do not indicate any negative evolution of the cluster mass, contrary to what is expected 
in a high density universe. This is consistent with results obtained by members of this group 
(Guerra & Daly 1996, 1998; Guerra, Daly, & Wan 1998; Daly, Guerra, & Wan 1998), who study 
the characteristic size of FRII sources and find the data strongly favor a low value for the density 
parameter O m , which is greater than 3a below unity. Note however, we only study the central 
regions of clusters. 



-19- 



The values of gas mass fraction obtained for the clusters surrounding the FRII sources 
studied here are consistent with observed values for central regions of clusters, after the correlation 
between the radio spectral index and redshift is taken into account, and a redshift evolution of 
the cluster core radius is considered. This suggests that the cluster core radius evolution and the 
effects of the radio spectral index-redshift correlation are important. 

The gas mass fraction seems to decrease slowly with redshift. The increase of cluster core 
radius with redshift can be one cause. Another possible explanation is that a large fraction of the 
gas is still in the galaxies at high redshift, and is later stripped from the galaxies into the cluster 
at low redshift. Cooling of the cluster gas also tends to cause the gas mass fraction in the cluster 
center to increase with time. 

The results on mass and gas fraction should be taken as preliminary because of the small size 
of the sample. 

It is a pleasure to thank Greg Wellman for helpful discussions. This work was supported in 
part by the US National Science Foundation, the Independent College Fund of New Jersey, and 
by a grant from W. M. Wheeler III. 
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Fig. 1. — The luminosity in directed kinetic energy of the jet vs. the core-hot spot separation. The 
best-fitted lines for all galaxies, all galaxies except Cygnus A, and all radio-loud quasars, are shown 
with a solid line, a dashed line and a dotted line, respectively. Slopes of the best-fitted lines are 
labeled. The number in parenthesis is the uncertainty on the slope times \/rx 2 , which includes the 
effect of a reduced x 2 > defined as rx 2 , greater than one for the fit. 

Fig. 2. — The luminosity in directed kinetic energy of the jet vs. redshift. 

Fig. 3. — The luminosity in directed kinetic energy of the jet vs. the radio power at 178 MHz. 

Fig. 4. — Lifetime of the collimated outflows vs. redshift. 

Fig. 5. — Thermal pressure of the ambient gas in the vicinity of the radio lobe vs. the core-hot 
spot separation. The best-fitted lines for all sources, including galaxies and quasars, and all sources 
except Cygnus A are shown with a dashed line and a dash-dot line, respectively. Only sources with 
detections on the Mach number M, and hence Pth, are included in the fit. 

Fig. 6. — Thermal pressure of the ambient gas in the vicinity of the radio lobe vs. the core-hot spot 
separation for two redshift bins (z < 0.9, z > 0.9). In the low-redshift bin, slopes of the best-fitted 
lines are -1.43 ± 0.21(0.27) with a reduced x 2 of 1.58 for all sources, and -1.35 ± 0.25(0.32) with 
a reduced x 2 of 1.65 for all sources except Cygnus A. In the high redshift bin, the slope of the 
best-fitted line is —0.41 ± 0.20(0.30) with a reduced x 2 of 2.29. Only sources with detections on 
Pth are included in the fit. 

Fig. 7. — Thermal pressure of the ambient gas at the cluster center vs. redshift for an isothermal 
gas distribution that follows the King density profile. A non-evolving cluster radius of 50/i -1 kpc, 
and Po = 2/3 are used. The same plots for an evolving core radius of 50(1 + z) 16 ^ 1 kpc are very 
similar to these and are not shown here. Bounds are not included in the fits. 

Fig. 8. — Thermal pressure of the ambient gas at the cluster center vs. the radio power at 178 
MHz for an isothermal gas distribution that follows the King density profile. A non-evolving 
cluster radius of 50/i -1 kpc, and 0o = 2/3 are used. The same plots for an evolving core radius of 
50(1 + z) 16 /i _1 kpc are very similar to these and are not shown here. Bounds are not included in 
the fits. 

Fig. 9. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas distribution 
that follows the King density profile with an evolving cluster core radius of 50 (1 + z) 1 - 6 /?,^ 1 kpc 
and (3q = 2/3 is used. The radio spectral index is redshift-corrected. Sources that are cooling 
at the point where the cluster temperature T is measured are marked with pentagons. The best- 
fitted lines for all sources, including both galaxies and quasars, all sources except Cygnus A, and 
all sources that are not cooling, are shown with a dashed line, a dot-dash line, and a solid line, 
respectively. Only sources with detections on T, hence mass, are included in the fit. 

Fig. 10. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas distribution 
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that follows the King density profile with an evolving cluster core radius of 50 (1 + z) ls h 1 kpc 
and fa = 2/3 is used. The radio spectral index is not redshift-corrected. 

Fig. 11. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas distribution 
that follows the King density profile with a fixed cluster core radius of 50 h" 1 kpc and (3q = 2/3 is 
used. The radio spectral index is redshift-corrected. 

Fig. 12. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas distribution 
that follows the King density profile with a fixed cluster core radius of 50/i _1 kpc and (3$ = 2/3 is 
used. The radio spectral index is not redshift-corrected. 

Fig. 13. — Gas fraction (M gas /M t ) within radius r vs. redshift. The radio spectral index is redshift- 
corrected, and an isothermal gas distribution that follows the King density profile with an evolving 
cluster core radius of 50 (1 + z) 1-6 /i -1 kpc and (3q = 2/3 is used. 
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Table 2. Results Obtained by Fitting M t (r) oc r ni (1 + z) n2 P\ 



type a 


No. b 


a-zT 


r c -zl A 




f 

ri2 


n 3 s 




N.C.C1. 


11 


Yes 


Yes 


2.38 ± 0.70(0.52) 


6.86 ± 5.48(4.1) 


-1.31 ± 1.57(1.17) 


0.56 


G+Q-C 


14 


Yes 


Yes 


2.95 ± 0.31(0.23) 


2.72 ± 3.80(2.84) 


-0.08 ± 1.01(0.76) 


0.56 


G+Q 


16 


Yes 


Yes 


2.88 ± 0.30(0.23) 


-0.73 ± 0.96(0.73) 


0.80 ± 0.37(0.28) 


0.58 


N.C.C1. 


13 


No 


Yes 


3.05 ± 0.36(0.28) 


4.46 ± 4.11(3.18) 


-0.13 ± 1.06(0.82) 


0.60 


G+Q-C 


14 


No 


Yes 


2.95 ± 0.31(0.23) 


3.66 ± 3.82(2.88) 


0.02 ± 1.02(0.77) 


0.57 


G+Q 


1G 


No 


Yes 


2.89 ± 0.30(0.23) 


0.49 ± 0.96(0.73) 


0.83 ± 0.37(0.28) 


0.58 


N.C.C1. 


11 


Yes 


No 


1.36 ± 0.63(0.52) 


7.40 ± 4.55(3.72) 


-0.93 ± 1.33(1.09) 


0.67 


G+Q-C 


14 


Yes 


No 


2.12 ± 0.27(0.23) 


3.18 ± 3.33(2.85) 


0.44 ± 0.90(0.77) 


0.73 


G+Q 


16 


Yes 


No 


2.11 ± 0.26(0.21) 


1.63 ± 0.91(0.75) 


0.84 ± 0.33(0.27) 


0.68 


N.C.C1. 


13 


No 


No 


2.17 ± 0.31(0.28) 


4.45 ± 3.40(3.10) 


0.47 ± 0.88(0.80) 


0.83 


G+Q-C 


14 


No 


No 


2.12 ± 0.27(0.23) 


4.09 ± 3.20(2.77) 


0.53 ± 0.86(0.74) 


0.75 


G+Q 


16 


No 


No 


2.11 ± 0.27(0.22) 


2.84 ± 0.91(0.76) 


0.86 ± 0.33(0.27) 


0.69 



a Type of sources included in the fit: "G+Q" refers to all sources with temperature estimates, including 
radio galaxies and radio-loud quasars, "G+Q-C" refers to sources other than Cygnus A, and "N.C.C1." refers 
to clusters that are not cooling at the position where the temperature is measured. 

b Number of data points used for the fit. 

c Whether the redshift-correction on the radio spectral index is applied. 

d Whether a redshift evolution of the cluster core radius r c is considered. The core radius is taken to be 
50 (1 + z) 16 /i _1 kpc if an evolution is considered, and 50 hr 1 kpc if no evolution is considered (see § |4.2| ). 

c ni and its error. The number in parenthesis is the error on ri\ times \J ry 2 which includes the effect of 
the reduced y 2 , defined as ry 2 , of the fit. 

' g same as note c for ni and ri3, respectively. 
h Rcduced y 2 of the fit. 
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Table 3. Results Obtained by Fitting M t (r) oc r" 1 (1 + z) n2 



type a 


No. b 


a-zl c 


r c -z? d 


rii c 


f 


9 

rx 


N.C.C1. 


11 


Yes 


Yes 


2.70 ± 0.59(0.45) 


2.39 ± 1.24(0.94) 


0.58 


G+Q-C 


14 


Yes 


Yes 


2.95 ± 0.31(0.22) 


2.43 ± 1.05(0.75) 


0.51 


G+Q 


16 


Yes 


Yes 


2.67 ± 0.29(0.28) 


0.63 ± 0.73(0.69) 


0.90 


N.C.C1. 


13 


No 


Yes 


3.04 ± 0.35(0.26) 


3.99 ± 1.16(0.86) 


0.55 


G+Q-C 


14 


No 


Yes 


2.95 ± 0.31(0.22) 


3.74 ± 1.06(0.76) 


0.52 


G+Q 


1G 


No 


Yes 


2.67 ± 0.29(0.28) 


1.90 ± 0.73(0.70) 


0.93 


N.C.C1. 


11 


Yes 


No 


1.60 ± 0.52(0.42) 


4.30 ± 1.04(0.84) 


0.65 


G+Q-C 


14 


Yes 


No 


2.13 ± 0.27(0.22) 


4.75 ± 0.83(0.68) 


0.68 


G+Q 


1G 


Yes 


No 


1.88 ± 0.25(0.27) 


3.34 ± 0.62(0.66) 


1.13 


N.C.C1. 


13 


No 


No 


2.20 ± 0.31(0.27) 


6.19 ± 0.94(0.82) 


0.77 


G+Q-C 


14 


No 


No 


2.13 ± 0.27(0.23) 


6.00 ± 0.83(0.70) 


0.72 


G+Q 


1G 


No 


No 


1.88 ± 0.25(0.27) 


4.57 ± 0.62(0.67) 


1.16 



same as those in Table 0. 
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10 



0.1 



solid: gal. 
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open: qua. 



i r 



i i i i i I i r 

sip: 0.61±0!09(0.29) x 2 - 9.16 

sip: 0.58±0. 1(0.3) xl- 9 - 64 

sip: 0.17±0. 13(0.67) x 2 - 27 - 38 
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TABLE 1 
Source Properties 



Source 




z 


LogL,(I) a 


LogL,(II) b 


LogP tA c 


LogP thc (A) d 


LogP thc (B) c 




oU 


r r 

55 


G 


0.72 


45.55 ±0.15 


45.39 ±0.15 


-11.34 ±0.19 


-10.26 ±0.19 


-10.87±0.19 


1. 










45.49 ±0.18 


45.33 ±0.18 


-11.08 ±0.19 


-9.98 ±0.19 


-10.60 ±0.19 


1. 


ou 


08.1 


Q 


1.238 


46.11 ±0.12 


45.88 ±0.12 


-10.50 ±0.19 


-9.61 ±0.19 


-10.32 ±0.19 


0. 










45.53 ±0.15 


45.28 ±0.14 


< -10.25 


< -9.26 


< -10.03 


1. 


3C 


68.2 


G 


1.575 


45.53 ±0.30 
45.73 ±0.17 


45.25 ±0.29 
45.48 ±0.17 








1. 
1. 


3C 


154 


Q 


0.58 


44.30 ±0.12 
44.65 ±0.12 


44.17 ±0.12 
44.51 ±0.12 








9. 
5. 


3C 


175 


Q 


0.768 


45.18 ±0.11 


45.02 ±0.11 


-10.50 ±0.15 


-9.78 ±0.15 


-10.27 ±0.15 


2. 










45.57 ±0.11 


45.41 ±0.11 


-10.85 ±0.09 


-9.91 ±0.09 


-10.50 ±0.09 


1. 


oU 


239 


G 


1.79 


46.16 ±0.13 


45.87 ±0.12 


-9.80 ±0.20 


-9.55 ±0.20 


-9.79 ±0.20 


0. 


oL 


247 


G 


0.749 


44.99 ±0.15 
44.69 ±0.16 


44.83 ±0.15 
44.54 ±0.16 


-10.54 ±0.57 


-10.45 ±0.57 


-10.53 ±0.57 


3. 
5. 


oU 


254 


Q 


0.734 


44.97 ±0.13 


44.80 ±0.12 


< -9.54 


< -9.33 


< -9.50 


3. 




265 


G 


0.811 


44.97 ±0.19 


44.80 ±0.19 


< -9.97 


< -8.79 


< -9.48 


3. 










45.26 ±0.19 


45.10 ±0.19 


< -9.86 


< -8.82 


< -9.47 


2. 


3C 


267 


G 


1.144 


45.58 ±0.19 


45.36 ±0.19 


< -9.16 


< -8.39 


< -9.01 


1. 










45.58 ±0.14 


45.38 ±0.14 


< -9.64 


< -8.92 


< -9.51 


1. 


3C 


268.1 


G 


0.974 


45.32 ±0.13 


45.18 ±0.13 


< -10.17 


< -9.30 


< -9.93 


2. 










45.66 ±0.12 


45.55 ±0.12 


< -10.00 


< -9.24 


< -9.81 


1. 


3C 


268.4 


Q 


1.4 


45.99 ±0.22 


45.72 ±0.21 


< -8.79 


< -8.62 


< -8.77 


1. 


3C 


270.1 


Q 


1.519 


45.96 ±0.19 


45.64 ±0.18 








1. 










45.96 ±0.14 


45.71 ±0.14 


< -8.76 


< -8.60 


< -8.75 


1. 


3C 


275.1 


Q 


0.557 


43.95 ±0.17 


43.80 ±0.17 


< -9.75 


< -9.53 


< -9.69 


15 










44.61 ±0.14 


44.47 ±0.14 


< -9.36 


< -9.25 


< -9.33 


6. 


3C 


280 


G 


0.996 


45.29 ±0.13 
44.98 ±0.11 


45.11 ±0.13 
44.79 ±0.11 








2. 
3. 


3C 


289 


G 


0.967 


45.15 ±0.11 


44.96 ±0.11 


-10.25 ±0.12 


-10.13 ±0.12 


-10.24 ±0.12 


3. 










45.18 ±0.12 


45.00 ±0.12 


< -9.35 


< -9.23 


< -9.33 


2. 



TABLE 1 Continued 



Source 




z 


LogL,(I) a 


LogL,(II) b 


LogP th c 

O III 




LogP thc (A) d 


LogP thc (B) e 


3C 322 


G 


1.681 


46.11 ±0.20 


45.79 ±0.19 








1 








46.14 ±0.12 


45.88 ±0.12 


-10.08 ±0. 


.14 


-9.47 ±0.14 


— 10.03 ± 0.14 i 


3C 330 


G 


0.549 


45.24 ±0.18 


45.16 ±0.18 


-11.19 ±0. 


.23 


-10.29 ±0.23 


-10.76 ±0.23 : 








45.15 ±0.15 


45.02 ±0.14 


-11.18 ±0. 


.28 


-10.30 ±0.28 


-10.76 ±0.28 : 


3C 334 


Q 


0.555 


A A C\ A 1 C\ "1 O 

44.94 ± 0.12 
45.05 ±0.22 


A A a i I n 1 n 

44.81 ± 0.12 
44.89 ± 0.20 










3C 356 


G 


1.079 


45.62 ±0.21 


45.41 ±0.21 


-10.51 ±0 


.22 


-9.14 ±0.22 


-10.01 ±0.22 








45.88 ±0.21 


45.69 ±0.21 


-10.59 ±0 


.17 


-9.54 ±0.17 


-10.29 ±0.17 


3C 405 


G 


0.056 


45.07 ±0.13 


45.06 ±0.13 


-10.31 ±0 


.19 


-10.07 ±0.19 


-10.10 ±0.19 ; 








44.99 ±0.12 


44.97 ±0.12 


-10.04 ±0 


.17 


-9.74 ±0.17 


-9.78 ±0.17 : 


3C 427.1 


G 


0.572 


44.72 ±0.13 


44.60 ±0.13 


< -9.81 




< -9.52 


< -9.72 








44.83 ±0.11 


44.71 ±0.11 


-10.37 ±0. 


.09 


-10.05 ±0.09 


-10.27 ±0.09 - 



a Logarithm of the luminosity in directed kinetic energy of the jet in h 2 erg s 1 . 

b Logarithm of the luminosity in directed kinetic energy of the jet in h~ 2 erg s _1 , with redshift-corrected 
(see §3.1). 

c Logarithm of the thermal pressure of the ambient gas in the vicinity of the radio lobe in /i 4 / 7 dyne < 

bounds are listed for source without dections. 

d Logarithm of the thermal pressure at the center of the surrounding cluster in /i 4 / 7 dyne cm" 2 , assuming 

of 50 k- 1 kpc (see §4.2). 

e Logarithm of the thermal pressure at the center of the surrounding cluster in h A l 7 dyne cm" 2 , assuming 
of 50 (1 + z) 1 - 6 h- 1 kpc (see §4.2). 

Characteristic lifetime of the outflow in (1 ±0.25) x 10 7 /i 2//3 yrs, where the uncertainty is that on the nc 
g Characteristic lifetime of the outflow in (1 ± 0.25) x 10 7 h 2 ^ 3 yrs, with redshift-corrected radio spectral 
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